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We consider the hydrodynamic theory of an active fiuid of self-propelled particles with nematic 
aligning interactions. This class of materials has polar symmetry at the microscopic level, but 
forms macrostates of nematic symmetry. We highlight three key features of the dynamics. First, as 
in polar active fluids, the control parameter for the order-disorder transition, namely the density, 
is dynamically convected by active currents, resulting in a generic, model independent dynamical 
self-regulation that destabilizes the uniform nematic state near the mean-field transition. Secondly, 
curvature driven currents render the system unstable deep in the nematic state, as found previously. 
Finally, and unique to self-propelled nematics, nematic order induces local polar order that in turn 
leads to the growth of density fluctuations. We propose this as a possible mechanism for the smectic 
order of polar clusters seen in numerical simulations. 



I. INTRODUCTION 

Active materials are soft materials driven out of equi- 
librium by energy input at the microscale. This liberates 
the fluctuations from the constraints of equilibrium such 
as fluctuation-dissipation relations and reciprocity. As 
a consequence, several exotic emergent behaviors result, 
such as long range order in 2D [U [2], anomalous fluc- 
tuations O E] , dynamical structures and patterns [MS] ■ 
In addition to serving as prototypical systems to explore 
emergent dynamical behavior, active materials also form 
the physical scaffold of biological systems in that active 
matter, when coupled to regulatory signaling pathways, 
provides a model for a variety of living systems, such as 
bacterial biofilms or the cytoskeleton of a cell. 

Active particles are generally elongated and form ori- 
entationally ordered states [5] . The nature of the ordered 
state depends on both the symmetry of the individual 
particles and the symmetry of the aligning interactions 
(see Fig. [T]). Physical realizations of polar active par- 
ticles (characterized by a head and a tail) include bac- 
teria, asymmetric vibrated granular rods, and polarized 
migrating cells. Polar active entities are often modeled as 
self-propelled (SP) particles, where the activity is incor- 
porated via a self-propulsion velocity of the individual 
entities. Apolar (head-tail symmetric) active particles, 
often referred to as " shakers" , have also been considered 
in the literature. Realizations are symmetric vibrated 
rods It has also been argued that melanocytes, the 
cell that distribute pigments in our skins, may effectively 
behave as "shakers" [101 E]. 

The nature of the interaction is of course crucial in con- 
trolling the properties of the ordered state. Apolar active 
particles generally experience apolar interactions and the 
resulting ordered state has the symmetry of equilibrium 
nematic liquid crystals. The broken orientational sym- 
metry identifies a direction n, but the ordered state is 
invariant for fi — >■ — fi. The properties of these active ne- 




FIG. 1. Top - active particles of various microscopic symme- 
try: (a) Polar active particles with head/tail asymmetry re- 
sulting in polar interactions, as studied in [13lll4j . (b) Apolar 
active particles, as studied in |3l [12] , (c) Sels-propelled par- 
ticles, with physical head-tail symmetry, resulting in apolar 
interactions. Bottom - ordered macroscopic states of active 
particles: polar active fluid (left) formed by polar particles (a) 
with polar interactions; active nematic fluid (center) formed 
by apolar particles (b) with apolar interactions; self-propelled 
active nematic fluid (right) formed by self-propelled particles 
(c) with apolar interactions. 



matic fluids have been studied by several authors [3l [12] . 
SP and polar particles may experience either polar in- 
teractions, i.e., ones that tend to align particles head to 
head and tail to tail, or interactions that are apolar, i.e, 
align particles regardless of their polarity. A well stud- 
ied example of polar particles with polar interactions is 
provided by Vicsek-type models [T31 [T3]. This class of 
active systems can order in polar states, characterized 
by a nonzero vector order parameter and mean motion, 
and will be referred to as active polar fluids. 
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One can envisage a third class of active fluids that con- 
sists of SP particles whose emergent macrostates are ne- 
matic in symmetry. A realization of this is self-propelled 
particles with physical interactions, such as steric repul- 
sion or hydrodynamic interactions among swimmers in 
a suspension. It has been shown that these interactions 
lead to nematic, rather than polar order [151115]. What 
is key is the fact that interactions such as steric colli- 
sions or hydrodynamic couplings individually conserve 
momentum and hence cannot lead to the development 
of a macroscopic momentum for the system. The or- 
dered state of these systems has nematic symmetry, but 
as we will see below its properties are distinct from those 
of active nematics composed of apolar particles. Here 
we refer to this third class of active systems as self- 
propelled nematic fluids. Further, it has recently become 
apparent that models of polar particles with apolar in- 
teractions may be relevant to a number of physical sys- 
tems, including gliding myxobacteria [T7], suspensions 
of self-catalythic Janus colloids [TB] and motile epithe- 
lial cell sheets such as those studied in wound healing 
assays [111 HO] . Self-propelled nematics therefore repre- 
sent an important new class of active systems of direct 
experimental relevance. 

A useful theoretical framework for describing the col- 
lective behavior of active systems is a continuum model 
that generalizes liquid crystal hydrodynamics to include 
new terms induced by activity [TH [?11E5] . In this pa- 
per, we consider a minimal phenonienological continuum 
model of self-propelled nematics based on the equations 
derived by us from a microscopic model of SP hard rods. 
We show that in these systems it is important to explic- 
itly retain the dynamics of both the collective velocity 
or polarization field of the particles and of the nematic 
order parameter to unfold the mechanisms at play in the 
formation of emergent structures. Further, we show that 
all active nematics (both consisting of shakers and SP 
particles) exhibit the phenomenon of dynamical self reg- 
ulation, due to the fact that the parameter controlling the 
order-disorder transition, namely the density p of active 
particles, is not externally tuned, as in systems undergo- 
ing equilibrium phase transitions, but it is dynamically 
coupled to the order parameter. This coupling is ana- 
logue to the one present in polar fluids [H (U [301 and is 
a generic mechanism for emergent structure in all active 
systems, as demonstrated in our recent work j31] . 

The layout of the paper is as follows. First, we con- 
struct the hydrodynamic description of a self-propelled 
nematic using symmetry considerations highlighting the 
key features that distinguish the dynamics of the system 
from that of an active nematic. Then we examine the lin- 
ear stability of the homogeneous nematic state. We show 
that there exists three dynamical mechanisms responsi- 
ble for emergent structures in active fluids with nematic 
symmetry. The first is a model-independent instability 
that occurs in the vicinity of the mean field order disorder 
transition due to the coupling between order parameter 
and mass transport which renders the dynamics of the 



system self regulating. We argue that this instability is 
the basis for the emergence of bands and phase separation 
found ubiquitously active systems [5^H5i] . The second 
is the well known instability of director fiuctuations that 
arises from noncquilibrium curvature-induced fiuxcs and 
is closely related to the giant number fluctuations ob- 
served in these systems [3J jT^] . These two instabilities 
are common to both active nematics and self-propelled 
nematics, i.e., occur regardless of the symmetry of the mi- 
crodynamics. Finally, we show that there exists a third 
instability unique to self-propelled nematic fluids due to 
the fact that in these systems, large scale nematic order 
can induce local polar order, which in turn destabilizes 
the density. This mechanism may be responsible for the 
smectic order of polar clusters observed recently in sim- 
ulations of SP rods [331 13S] . We conclude with a brief 
discussion. 



II. THE MACROSCOPIC THEORY 

The hydrodynamic equations of a self-propelled ne- 
matic have been derived from systematic coarse-graining 
of specific microscopic models [TS1|T5]. Here we introduce 
these equations phenomenologically, with the goal of ex- 
amining the dynamics without the limitations imposed 
by the specific parameter values obtained from a micro- 
scopic model or resulting from the choice of the closure 
used in the kinetic equation. 

We limit ourselves to overdamped systems in two di- 
mensions. The hydrodynamic equations are then writ- 
ten in terms of three continuum fields: the conserved 
number density p (r, t) of active units, the polarization 
density T{r,t) — p(r, t)P (r, i), with P(r,i) a polariza- 
tion order parameter, and the nematic alignment density 
tensor Qij{r,t) = p{r,t)Sij (r,i). The polarization P is 
directly proportional to the collective velocity of the ac- 
tive particles, while Sij is the conventional nematic order 
parameter tensor familiar from liquid crystal physics. For 
a uniaxial system in two dimensions, Qij is a symmet- 
ric traceless tensor with only two independent compo- 
nents and can be written in terms of a unit vector n as 
Qij — Qifiifij — -^Sij), where Q = pS; S is the magnitude 
of the order parameter and the director n identifies the 
direction of spontaneously broken symmetry in the ne- 
matic state. For simplicity most of the discussion below 
refers to the case where the active particles are modeled 
as long thin rods with repulsive interactions. 



A. Active Nematic Hydrodynamics 

We first construct the dynamical equations of the sys- 
tem including single-particle convection terms induced by 
self propulsion, but assuming that self-propulsion does 
not modify the interaction between two rods. The hy- 



drodynamic equations then take the form 
dtp + voV ■ T = DV'^p 

dtT + Dr-T + UqV • Q 
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(la) 
(lb) 



^tQ^J-Dr [a (p) - /3Q:Q] g„- + v^F,, = DbV'g, 



Dsdk (diQkj + djQi,k - SijdiQ 



kl) 



(Ic) 



where Q:Q — QkiQki, is the rotational diffusion rate, 
and Fij = (diTj + djTi — SijV ■ r). All terms propor- 
tional to vq arise from one-particle convection due to 
self-propulsion and are the only consequence of activity 
in this simple model. The repulsive interactions among 
the particles generate the cubic homogeneous term (with 
/3 > 0) on the right hand side of Eq. (Ic) and a change 
in sign of a{p) ^ p — pc at a critical density pc, control- 
ling the transition between the isotropic and the nematic 
states. ^37j Interactions also give density-dependent cor- 
rections to the various diffusion coefficients for density 
(D), polarization (Dr), splay (Ds) and bend (Df,) defor- 
mations of the nematic alignment tensor. We will ignore 
all such corrections in the following [3H|- Due to the fact 
that the interactions are purely nematic, the polarization 
T decays on short time scales ^ D^^ for all strengths of 
activity. At long time, a hydrodynamic description can 
then be obtained by neglecting dtT in Eq. (lb), and using 



Eq. ( lb I to eliminate r from the other equations, with 



the result (to leading order in gradients) 
dtp = DV^P + VqVV:Q 

dtQ^j-Dr [a [p) - /3Q:Q] Q,j - /^^V^ 



(2a) 



Dsdk [diQkj 
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djQik - SijdiQki) 



+ Vp{d,d,-^-5,j^^)p 



(2b) 



This procedure yields qualitative new terms proportional 
to "Dq and "Dp that vanish in equilibrium. The term pro- 
portional to I?Q is the curvature induced density flux 
that has been discussed extensively by Ramaswamy and 
collaborators [3] and shown to be responsible for giant 
number fluctuations in the ordered state of active ne- 
matic. The diffusive coupling proportional to 'Dp de- 
scribes similar physics but has not been considered in 
earlier description of active nematic fluids. In addition, 
activity also yields corrections to the various diffusion co- 
efficients. We have, however, implicitly neglected those 
here by retaining the same notation for these quantities 
as in Eqs. (lal-(lc) to highlight the difference between 
these corrections that do not change the dynamics qual- 
itatively and the new terms proportional to TDq and Dp. 
Although obtained here by considering a system of self- 
propelled particles, Eqs. (2a) and (2b) have the same 



structure as the hydrodynamic equations of an active ne- 
matic, consisting of a collection of apolar active particles 
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FIG. 2. Illustration of momentum conserving colli- 

sions among self propelled particles. It can readily be 
shown that two rods as shown in (a), coming in with 
only their self-replenishing velocities, will acquire oppo- 
site angular momenta oji ~ z£vo [z ■ [u]-^ x U2) and u}2 ~ 
—zivo (ui • U2) [z • (ui X U2)], where the vectors are defined 
in the image and in |40|. The collision will therefore induce 
rotations as indicated, promoting alignment of the two rods. 
On the other hand, two nearly antialigned rods as in (b) ac- 
quire angular momenta of the same sign, inducing rotation of 
both rods in the same directions, and leaving their relative 
angle unchanged. 



(shakers) with apolar interactions. This is an important 
point as it stresses that the qualitative differences be- 
tween active and self-propelled nematic that have been 
observed in simulations must arise entirely from the de- 
pendence of the interaction on self propulsion vq. 



B. Momentum-conserving interaction of 
self-propelled nematogens 

As shown in Ref. |15j and supported by simulations 
of self-propelled hard rods [33 EH self-propulsion 
does modify the repulsive interaction in a qualitative way. 
This modification results in local build-up of polariza- 
tion in the nematic state, making it necessary to retain 
the dynamics of polarization density in the continuum 
model. The modification of the Onsager excluded vol- 
ume interactions among hard rods due to self propulsion 
is worked out in Ref. [40^. Here we simply give a qual- 
itative description of this effect and we refer the reader 
to that work for the technical details. First we note that 
the presence of a self propulsion speed along the long 
axis of the nematogen, results in a breaking of the ne- 
matic symmetry of the the microdynamics, as shown in 
Fig. [1] On the other hand, since the interactions conserve 
momentum, this cannot lead to a macroscopic breaking 
of polar symmetry as this would amount to the appear- 
ance of a spontaneous macroscopic momentum from a 
zero momentum state. Hence, only a homogeneous or- 
dered nematic state can occur and the associated mean 
field transition will be the same as in the case of the ac- 
tive nematic considered above, albeit with coefficients a 
and /3 renormalized by self- propulsion ^15, . Even though 
the polar symmetry cannot be broken macroscopically, 
momentum conservation allows the nematic ordering to 
induce local polar order in the system. To illustrate this, 
let us consider hard rods in two dimensions undergoing 
energy-momentum conserving interactions. As shown in 
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Fig. ([2]), the angular momentum transfer due to the lin- 
ear momentum from self-propulsion for a collision be- 
tween two rods scales as w ^ cos {9i — O2) sin {9i — 02)- 
If the rods arc nearly aligned head to head (as in Fig[2ja), 
the effect of this angular momentum is to turn the rods 
towards each other, while if they are nearly aligned head 
to tail as in Fig[2jb, the collision turns both rods in the 
same direction, leaving their relative angle unchanged. 
This mechanism effectively promotes head-to-head align- 
ment. Since collisions among such nearly aligned nemato- 
gens will dominate the dynamics in the nematic state, the 
nematic order effectively induces polar order. 



C. Self-propelled Nematic Hydrodynamics 



The fact that interactions among self-propelled ne- 
matogens tend to induce polar order is reflected in the 
hydrodynamic description by a number of new nonlinear 
terms that couple t and Q, with coefficients that vanish 
in the limit uq = 0. The continuum equations for a self- 
propelled nematic that incorporate the above physics are 
given by 



dtp + woV • T = D^'^p + PqVV : Q (3a) 
dtT + DrT + 7iQ:Qt - • Q + Ait • Vt= - i^qV • Q - y Vp + AarV • r + y Vr^ -I- D^VV (3b) 

^tQ^3 - Dr{a - /3Q:Q)Q»j + v^F,, + = Vp{d,dj - ^V^)^ + D^dk {d.Qkj + djQ^k ~ kAQki) + AV^Q^ 

(3c) 



where again we have implicitly neglected active correc- 
tions to D, Dt, Dg and Z?;, to highlight the new, purely 
active terms. Activity enters in Eqs. ([3| through the 
convective terms proportional to uq, the new terms with 
coefficients 7^ and to Ai, which vanish in equilibrium, 
as well as the terms proportional to TDq and 'Dp that 
arise here from active corrections to interactions. Fi- 
nally, the parameters a and j3 controlling the mean field 
isotropic-nematic transition are also renormalized by ac- 
tivity. The homogeneous nonlinearities proportional to 
7i in the polarization equation encode the fact that ne- 
matic order induces polar order. The latter is, however, 
only local as the equations do not a admit a homogeneous 



solution with nonzero r. Further, the active modifica- 
tion of the interactions, yield the convective nonlineari- 
ties ~ 0{tVt) that play a central role in the emergent 
physics of active polar fluids. 

Since the goal of this presentation is to highlight the 
mechanisms responsible for emergent structures, we sim- 
plify the equations by setting all of the equilibrium-like 
diffusion coefficients to be equal, i.e., D — — D\, — 
Dq, with the exception of the splay relaxation constants 
Dg. In addition, we assume A^ = A for all I's and 
7i = 72 = 7- Finally, we measure time in units of l/I?,- 
and lengths in units of The hydrodynamic 

equations then simplify to (in nondimcnsional form) 



dtp + vV -T = \l'^p + DqW : Q 
dtT + (1 + 7Q:Q) T - 7T • Q + At • Vr = -vV ■ Q - ^ Vp + A ( rV • x + ^\Jt^ ) + V^r 



(4a) 
(4b) 



dtQ^, -{a- /3Q:Q) Q,^ + vF,, + AG.j = Dp{d,d, ~ -5,jV^)p + Dsdk {d,Qk, + d^Q.k - S^.diQu) + V^Q,^ (4c) 



with Dq = Vq/Dq, Dp = Vp/Do and Ds = Ds/Dq. 
Finally, we assume a = — — 1 and /3 independent of 

Pc 

p. The effect of activity is assumed to affect the mean 
field phase transition only through the dependence of the 
critical density pc on the magnitude of self-propulsion 
speed. In this simplified form, the dynamics of the system 
is characterized by two central parameters: the mean 



density po of active nematogens and the self-propulsion 
velocity v = Vq/ VD^^Dqi which is effectively the Peclet 
number for this flow. The other parameters 7, A, Dq, 
Dp and Dg are in general functions of po and v, although 
we will treat them here as independent parameters and 
simply fix their values. 
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III. LINEAR DYNAMICS AND EMERGENT 
STRUCTURES 

The dynamics of self propelled rod-like particles with 
steric repulsion has been studied extensively by numeri- 
cal simulation of microscopic models |321|331[3SI[1T]. This 
work has revealed a rich variety of emergent structures, 
including bands of high density regions where the par- 
ticles are ordered along the direction of the band, lane 
formation, migrating defect structures and low Reynolds 
number turbulence. Here we examine the minimal con- 
tinuum model of self-propelled nematic given by Eqs. Q 
to identify the generic dynamical mechanisms responsi- 
ble for the emergence of these structures. As mentioned 
above, there are three important mechanisms for dynam- 
ical instabilities and parttern formation in these systems. 
To unfold the role of each of these mechanisms in control- 
ling the large-scale dynamics of the system, we analyze 
the linear stability of the ordered nematic state in various 
special cases that best highlight a particular mechanism. 

The ordered nematic state has constant density po, 
zero mean polarization density, Tq = 0, and a finite value 
for the nematic alignment tensor. Choosing a coordinate 
system with the x axis pointing along the direction of bro- 
ken nematic symmetry, the alignment tensor in the uni- 
form nematic state has components Qjjj. = —Qyy = Qo/2 

and Q°y = = 0, with Qo = \/ao//3 and oq = a{po). 
We now examine the linear stability of this state in var- 
ious regions of parameters by considering the dynamics 
of small fluctuations, Sp{r,t) = p{r,t) — po, ST{r,t) = 
T{r,t) and 6Qij{r,t) = Qij{r,t) — Qij- We will generally 
work in Fourier space by introducing Fourier transforms 
of the fluctuations as (p^it) = e^^''^S4>a(j,t), where 
(5^Q = {Sp,T,SQij). 



A. Dynamical Self-Regulation and Banding 
Instability 

We first consider the linear dynamics of the system in 
the region just above the mean-field transition at pc- For 
simplicity we only discuss spatial variations normal to 
the direction of broken symmetry, as these correspond to 
the most unstable modes, i.e., let k = fcy. Fluctuations 
in Tx and SQxy then decouple and are always stable. The 
dynamics of fiuctuations in Sp, Ty and SQyy is governed 
by three coupled equations. Fluctuations in Ty are always 
quickly damped near the mean-field transition, while the 
decay rate of 6Qyy, controlled to leading order by ao, 
vanishes as po ~^ pt ■ We therefore neglect fiuctuations 
in Ty and simply examine the coupled dynamics oi 5p and 
SQ = SQyy, given by 



dtSpk = -k'^Spk - Dqk^SQ]^ 
1 



(5a) 



2 



{l + Ds)k 




FIG. 3. The banding instability that occurs due to the 
self-regulating nature of the flow. The striped region is the 
parameter space in which this instability occurs. 



where a' — 



( §^ ) , or a' = 1/pc with the chosen pa- 

rameters. The decay of density and ordered parameter 
fiuctuations is then controlled by two coupled hydrody- 
namic modes. One of the modes has a finite decay rate 
(proportional to ao) in the limit fc — > and is always 
stable. At small wavevector, the dispersion relation of 
the other mode is given by 



(6) 



with S2 = 1 — -7^ 
where ao ^ 0, S2 < and S4 ^ 



(a'Qo + Dpk^)5pv 
(5b) 



and S4 > 0. Near the transition 

2 '2 

As a result, 

Sy{k) > for a range of wavevectors, resulting in the 
unstable growth of density and order parameter fiuctua- 
tions illustrated in Fig. [3 } Th e fastest growing mode has 
wave vector ko — S2/2S4 ^ (po — Pc)'^^^- Including 
the coupling to Ty will yield finite Peclet number cor- 
rections to the instability. Note that this instability is 
strongest in the vicinity of the order-disorder transition 
and is a manifestation of the fact that the dynamics of 
the system is self regulating, i.e., the control parameter 
associated with the phase transition, namely the density 
is dynamically coupled to the emergent ordering that re- 
sults from the transition through the curvature induced 
fiuxcs. This is the dynamics that leads the system to be 
intrinsically phase separated 

We recall that polar active fiuids exhibit a similar in- 
stability for wave vectors parallel to the direction of mean 
order. In that case the mode that goes unstable is a 
propagating mode and the instability signals the onset 
of solitary waves consisting of alternating ordered and 
disordered bands extending in the direction normal to 
that of mean order and traveling along the direction of 
broken symmetry. These bands have been observed in 
simulations of the Vicsek model [131 133 j as well as in 
numerical solutions of the nonlinear hydrodynamic equa- 
tions for polar fluids [43^-^45^. We have shown here that 
active nematics exhibit a similar instability, controlled 
by the interplay of of curvature currents (Dq) and the 
self- regulation due to the density dependence of a. The 
instability occurs even for v = 0, i.e., is present in both 
active and self-propelled nematics. It occurs for wavevec- 
tors perpendicular to the direction of broken nematic 
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symmetry and the mode that goes unstable is a diffu- 
sive one. It is therefore tempting to associate it with 
the emergence of the stationary bands consisting of al- 
ternating ordered (nematic) and disorders regions that 
have been seems in simulations of active systems with 
apolar interactions j32j and physical excluded volume 
interactions |33 ) .|35 ) . Finally, but most importantly, this 
instability mechanism is generic, in the sense that it does 
not depend on microscopic parameters, but only on the 
presence of a dynamical feedback between density and 
active currents. 



B. Curvature Induced Flux 

Next we consider the region of small w, A and 7, deep 
in the nematic phase. In this case, the long-wavelength 
dynamics is controlled by hydrodynamic modes associ- 
ated with fluctuations in the density and the director n. 
This case has been considered in the literature already 
and is summarized here for completeness O |36l l46] . 
For our choice of coordinates to linear order we have 



5Qa 



'5Qy 



and 



SQyx = QoSh{r,t). 



Neglecting polarization fluctuations that decay on mi- 
croscopic time scales, the linearized equations are given 
by 

dtSpw — —k'^Spk — QoDk^ sin26'(5fik , (7a) 

Dp — 
dtSfi^ = sin26'(5pk + [Ds + cos 261] PSfii, , (7b) 

where 6 is the angle between k and the direction of broken 
symmetry (x). 11 9 = 0,7r, the two equations are decou- 
pled and the modes are diffusive and stable. For gen- 
eral 9 one of the hydrodynamic modes becomes unstable 
for lyDpSin^ 29 > 2(Ds + cos 26*). This can be satisfied 
provided DqDp > 2Dg, i.e., the curvature driven fluxes 
exceed the restoring effects of diffusion. This instability 
has been discussed in detail elsewhere l36l. 



C. Induced Polar Order 

Finally, we examine the effect of fluctuations with spa- 
tial variations along the direction of broken symmetry, 
i.e., k = fcx. The relevant coupled fluctuations in this 
case are 5p, and 5Qxx- For simplicity, we consider the 
regime of large Peclet number iJ, where the linear dynam- 
ics is controlled by Euler order terms and neglect terms 
quadratic in the gradients, with the result 



(8) 



dtTxM + leTx,\L = -ikvSQxxM ~ ik-Spk (9) 

dtSQxx,k - '^SQxx,k = ik 2^0^ '^^'^ + ~^QoSpk 

(10) 



where 7e = 1 + 2*3o ~ i^o- As discussed earlier and 
highlighted in Fig. [2j the anisotropy of small angle colli- 
sions in the nematic state enhances polar order by sup- 
pressing the decay rate of Tx^k from its bare value of 1 
(in units of D~^) to 7e. The dispersion relations of the 
hydrodynamic modes associated with Eqs. ([6| are easily 
calculated at small wavevectors. Clearly, if 7e < 0, In 
addition, as a consequence of this built up of polar order, 
the diffusive mode associated with the density fluctua- 
tions (i.e., the only truly hydrodynamic mode in these 
considerations), given by 



Sx{k) 




(11) 



becomes unstable. In general 7e depends on microscopic 
details of the model, but there is no reason to exclude 
a priori that it could change sign and indeed does for 
the case of long thin hard rods with excluded volume 
interactions [36j . The linear analysis here is of limited 
utility because of the existence of a homogeneous insta- 
bility but is shown here to indicate that the build up of 
polarization due to the momentum conserving nature of 
the interactions has a dramatic consequence on the dy- 
namics of the system. This may indeed be the mechanism 
responsible for the smectic order observed within a single 
polar cluster in simulations of self-propelled rods [331 ES] ■ 
Finally, we stress that the nonlinear homogeneous terms 
proportional to 7 and responsible for the renormalization 
of 7e always vanish in an equilibrium state because the 
nematic symmetry of such a state by definition forbids a 
nonzero uniform value of the mean polarization. 



IV. DISCUSSION 

We have considered in this paper the hydrodynamics 
of active overdamped fluids that can order in nematic 
states. These are collections of active particles that in- 
teract via apolar (nematic) aligning interactions, such as 
steric repulsion or medium-mediated hydrodynamic cou- 
plings. One can identify two classes of such fluids, de- 
pending on the properties of the individual active units. 
Active nematics consist of shaker particles that are them- 
selves apolar. Self-propelled nematics are collections of 
particles that are physically head-tail symmetric (such 
as SP rods) , but where a microscopic dynamical polarity 
is induced by self-propulsion . Although both systems 
form ordered states of nematic symmetry, their dynam- 
ical behavior is qualitatively different, as seen in recent 
simulations 1331151 155]. 

The hydrodynamic equations of active nematics have 
the form given in Eqs. 2. We have shown that the same 
equations are also obtained by considering SP particles 
and neglecting the effect of self-propulsion on the inter- 
action between active units, suggesting that the active 
nematic may be considered the zero Peclet number v 
limit of self-propelled nematic. In this case the only ac- 
tive term is the curvature current proportional to Vq 
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in Eq. ( |2a[ ). This non equilibrium coupling of orienta- 
tion and flow induces instabilities of the ordered state 
that have been studied before in the literature [3 [TH [35] 
and are also summarized in section. IIII Bl The curvature 
current is also key in controlling the banding instability 
arising from dynamical self- regulation discussed in sec- 
tion |TTrX] In fact this instability, although not discussed 
before in the literature for overdamped active nematic, 
occurs in all active fluids of nematic symmetry, both for 
shakers and self-propelled particles. It arises from the 
density dependence of the parameter a{rp) that controls 
the mean-field transition and the fact that in active sys- 
tems p is not tuned from the outside, as in equilibrium, 
but is itself a dynamical variable convected by the order 
parameter. 

The hydrodynamic equations of self-propelled nemat- 
ics given in Eqs. ^ (or Eqs. Q in the dimensionless form 
studied here) contain many new active terms that arise 
from modifications of the two-body interaction due to 
self propulsion. These equations have also been derived 
by us for a specific microscopic model of self-propelled 
hard rods [TH |3D] , although the low order closure of the 
kinetic theory used in that work only gives terms up to 
quadratic in the hydrodynamic fields. Self-propelled ne- 
matics also exhibit both the curvature induced instability 
discussed in section p^IIB| and the banding instability dis- 



cussed in section UlI Al Both are of course modified at fi- 
nite Peclet number due to additional convective contribu- 
tions to the underlying mechanisms the details of which 
will be discussed elsewhere. In addition, self-propulsion 
yields a novel instability due to the built-up of local po- 
lar order discussed in section lTlI CI This arises because in 
the nematic state most binary collisions involve nemato- 
gens that are nearly aligned or anti aligned, as shown 
in Fig. [2] When the nematogens are self-propelled, col- 
lisions of nearly aligned and nearly anti aligned pairs 
are not identical. Nearly aligned pairs tend to further 
align upon collisions, while nearly anti-aligned pairs are 
turned away from each other. As a result, local polar or- 
der is enhanced and the nematic state becomes unstable 
as discussed in section |III C| It is tempting to associate 
this instability with the onset of "polar clusters" that 
have been observed ubiquitously in simulations of self- 
propelled rods [33l EH |39] , as well as in experiments in 
gliding myxobacteria [T7] . 
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